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Interesting magnetic properties and spin-exchange interactions along various possible pathways 
in half-integral spin quantum magnetic tetramer system: A2PO4OH (A = Co, Cu) are investi- 
gated. Interplay of structural distortion and the magnetic properties with the evolution of localized 
band structure explain the gradual transition from a three-dimensional antiferromagnet to a low- 
dimensional frustrated magnetic system along the series. Detailed study of exchange mechanism 
in this system explores various possibilities of complex magnetic interaction. Electronic structure 
of this series, studied with the help of different appropriate density-functional approaches like N-th 
order Muffin-Tin orbital(NMTO) and Planewave Pseudopotential calculations incorporating onsite 
Coulomb repulsion(U), identifies the underlying magnetic exchange mechanism of this series. There- 
after a generalized minimal model spin-Hamiltonian is constructed for the low-dimensional system. 
Solution of this model Hamiltonian within first order perturbation theory results into the evaluation 
of spin-gap in the spin-tetramer system. In addition, the effects of size-confinement and volume 
reduction on the relevant exchange integrals and spin-gap of the low-dimensional system are also 
discussed. 

PACS numbers: 75.50.Ee, 71. 15. Mb, 75.10.Jm 



I. INTRODUCTION 

Thorough understanding of the interrelated structural 
and magnetic properties of complex geometrically frus- 
trated systems remain always interesting as absence of 
long-range magnetic order in such systems opens up the 
possibilities of a wide range of excited state properties 
PQ. Such half-integral spin, especially the spin- 1/2 low- 
dimensional magnetic system with a singlet ground state, 
attracts much attention due to their fundamental quan- 
tum nature [2] [3] and are also considered to be suitable 
two-level system for quantum computing. [4] Observation 
of pseudo spin-gap in high-T c cuprates [5] has motivated 
intensive study of spin-exchange mechanism of new spin- 
gap systems. Increasing experimental activity on com- 
plex low-dimensional magnetic systems 6 9 motivates 
a detailed theoretical first-principles investigations of the 
underlying mechanism of substitution-induced structural 
distortions and the resulting changes in magnetic ex- 
change properties to understand the evolution of a long- 
range ordered system into a low-dimensional one. More- 
over, for complex geometrically frustrated system, iden- 
tification of exchange paths relevant to the formation 
of minimal spin-Hamiltonian and thus determination of 
effective magnetic mechanism can provide the insight 
about magnetostructural correlations. 

A2PO4OH (A = Co, Cu), the natural mineral Li- 
bethenite, possessing emerald green and dark purple 
colour for Cu and Co based systems respectively, belongs 
to the phosphate mineral class p~0] and is mostly used 
for their widespread catalytic properties. [TT] The par- 
ent compound, C02PO4OH is experimentally observed 
to possess a three-dimensional antifcrromagnetic order 
at around 70 K. Substitution of Cu gradually in place of 
Col, Co2 and finally for both Col and Co2 sites leads to 



the onset of low-dimensional magnetic properties in this 
system. [HI [12] The final compound CU2PO4OH, a spin- 
l/2(<i 9 ) low-dimensional tetrameric system, being the 
most interesting among such class, will be investigated 
in more detail in the present work. Previous experimen- 
tal observations on CU2PO4OH system [T31 Q3] predict 
a singlet ground state with a spin-gap ^141 K. A sim- 
pler analogue of the spin-exchange mechanism were put 
forward via the square spin-tetramer(SQST) model. [14] 
In the present study, we adopt suitable Density Func- 
tional Theory (DFT)-based methodologies to first iden- 
tify and then calculate the relevant exchange interactions 
in such systems, which appears to be more complicated 
than the SQST model. [2] Next, we construct the mini- 
mal Heisenberg Model for the final system and approxi- 
mately solve it within first order perturbation theory to 
evaluate the spin-gap between the singlet ground state 
and the triplet excited state. Such multistep analysis 
will be helpful to build a generalised spin-model Hamil- 
tonian for the low-dimensional tetramer systems and also 
provides understanding of the correlation between the 
magnetic and structural properties. For the last system, 
effects of size-reduction and cell-volume reduction on the 
exchange mechanism are also studied. 



II. FIRST PRINCIPLES ANALYSIS 

A2PO4OH belongs to the space group Pnnm(58) with 
two types of the transition metal cation A (Al or A2, 
suitably replaced by Cu or Co) and four types of O- 
ligand(01,02,03 and 04). Al forms a distorted octa- 
hedra with 01, 02 and 04(H), whereas A2 is fivefold 
coordinated with 01, 03 and 04(H)(Fig 1 and 5(a)). 
The experimental value of lattice parameters are ob- 
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FIG. 1: (colour online)The colour code of the atoms are 
Cu(Co) - Cyan, O - red, P - magenta, H - grey. (a): Corner 
shared A106 octahedra and A2O5 trigonal bipyramids are 
shown for comparison with the actual structure, (b): Actual 
atomic coordination showing A104(OH)2 edge-sharing infi- 
nite octahedral chains along the z-axis and A204(OH) edge- 
sharing trigonal bi-pyramidal dimers. 



tained from reference |10j . For a detailed understand- 
ing of magnetic interactions in the tetrameric system, we 
have investigated all four systems in the series, viz. (I) 
C02PO4OH, (II) CuCoP0 4 OH(Cu replacing Col), (III) 
CoCuP0 4 OH (Cu replacing Co2) and (IV) Cu 2 P0 4 OH, 
all generated from the parent Libethenite system I. First 
principles Density functional investigation of this series 
involves three steps; viz. (1) For all four systems, lattice 
parameters and the atomic positions in the 36-atom unit 
cell are optimized using conjugate-gradient relaxation 
criteria as implemented in VASP [T5J [TB] with projector 
augmented wave (PAW) [17] formalism. Bond-length and 
bond-angles of the relaxed structure provide a prior idea 
of the effects of structural distortion on the electronic 
and magnetic properties due to Cu substitution; (2) With 
the optimized positional coordinates and lattice parame- 
ters as input, we have utilized N-th order Muffin-Tin or- 
bital (NMTO) downfolding technique to obtain the tight- 
binding hopping parameters for the Cu and Co localized 
rf-bands. Under second-order perturbation theory, com- 
putation of these hopping paramaters (ii)corresponding 
to some specific paths may lead to a quantitative com- 
parison of the exchange couplings along them. For an- 
tiferromagnetic couplings, 4f?/J7 is a measure of the ex- 
change coupling Ji, with U being the onsite Coulomb 
interaction strength for Co and Cu - 3c? bands; (3) After 
obtaining these relevant exchange paths, we have esti- 
mated the exchange coupling constants and the nature of 
exchange along these pathways by a total energy calcula- 
tion for different spin-configurations in various supercell 
configurations consisting of 72 atoms within Local Spin- 
density Approximation(LSDA) + onsite Coulomb repul- 
sion^) formalism with VASP-PAW potentials. The first 
step consists of the analysis of optimized structure. With 
gradual substitution of Cu in the parent C02PO4OH sys- 
tem, the relaxed structure indicates that due to smaller 
size of Cu 2+ ion, there is sizable amount of structural 
distortions including change in the bond angle and bond- 
length around the site of substitution. The trend of 



TABLE I: Important intra and inter-tetramer distances (d) 
and transition-metal cation to ligand bond-lengths (b)and 
bond-angles for the relaxed structures of the Libethenite- 
series. See fig 5(a) for structure. 



d(A) 


I 


II 


III 


IV 


A1-A2 (intra) 


3.47 


3.43 


3.35 


3.30 


Al-Al (intra) 


2.78 


2.82 


2.90 


2.95 


Al-Al (inter) 


2.89 


2.91 


2.89 


2.95 


b(A) 


I 


II 


III 


IV 


Al-Ol 


1.91 


1.93 


2.00 


1.97 


A2-01 


2.03 


1.97 


1.87 


1.93 


A2-03 


1.95 


1.93 


1.99 


2.14 


Al-04 


2.28 


2.40 


1.98 


2.32 


Angle (degrees) 


I 


II 


III 


IV 


02-A1-02 


96.52 


90.35 


86.88 


82.47 


Ol-Al-Ol 


83.88 


83.03 


82.50 


81.20 


A1-02-A1 


95.99 


96.64 


97.99 


98.40 


Al-Ol-Al 


96.11 


96.33 


98.09 


97.35 


A2-03-A2 


99.71 


100.78 


102.82 


101.42 


A1-Q1-A2 


123.40 


122.86 


122.30 


121.64 



significant bond-length, ionic distance and bond-angle 
modifications are presented in table I for comparison. 
From system I to IV, as a combined effect of decrease 
in 02-A1-02 and Ol-Al-Ol and increase in Al-02-Al 
bond-angles, the difference between the intra and inter 
tetramer Al-Al z-direction distance nullifies along the 
series (Table I) and thereby the z-direction Al-Al dimer 
chain structurally turn into a monomer chain. This also 
leads to a Jahn- Teller like elongation of the distorted oc- 
tahedra. Although the A1-01-A2 angle does not change 
much, but as a result of increase in Al-Ol-Al and A2- 
03- A2 angle, the intra-tetramer A1-A2 distance reduces 
and inter-tetramer A2-A2 distance increases along the 
series. Interestingly, disparity of bond-lengths of Al and 
A2 with 01 and 02 and H-ligands leads to a competition 
in the cationic covalency which also has footprints on the 
ligand-ficld splitting and localized band structure as dis- 
cussed in the next step. However, system (III) slightly 
breaks the gradual order in some structural details, viz, 
Al-Ol-Al and A2-03-A2 bond angles are more than sys- 
tem IV, which leads to a change in the coordination of 
ligand 04 around Al (dist ~ 1.98A) such that Al - 04 
bonds are formed for system III. 

Incorporation of H-ion into the ligand coordination sig- 
nificantly changes the ligand field and the resulting split- 
ting for this series of materials. The three-dimensional 
real structure of the system consists of A104(OH) 2 plane- 
sharing infinite octahedral chains along the z-axis, edge- 
shared with A204(OH) trigonal bi-pyramidal dimers. 
This structure is distorted from the edge-shared AlOg- 
octahedra, corner-shared with A2O5 trigonal bipyramids 
as shown in fig 1(a). Figures 1(a) and (b) depict a com- 
parison of the distorted ligand coordination around Al 
and A2 due to H-ion. The ligand field splitting for Al 
and A2 3d-levels are also deviated from the usual octa- 
hedral and trigonal-bipyramidal ones and becomes more 
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FIG. 2: (color online) The figure represents the LSDA band structure for all four systems (a)Co2P040H, (b)CuCoP040H(Cu 
replacing Col), (c) CoCuP0 4 OH (Cu replacing Co2) and (d) Cu 2 P0 4 OH with the localized NMTO-d-bands near Fermi-level 
in blue. Band-width reduction of the localized manifold from system I and II to III and IV can be observed. 



closer to a distorted square-antiprismatic (coordination 
no. 8) and octahedral (coordination no. 6)coordination 
for Al and A2, respectively. A closer look to fig 1(b) also 
reveal that actual coordination around A2 is distorted 
octahedral and two plane shared octahedra around Al 
along z resembles more of a distorted antiprism. This 
description will be more evident in the orbital-analysis 
of the different 3d-characters. 

In the second step, the relaxed structural inputs are 
used to extract the resulting effects of structural distor- 
tion on the localized bandstructure. Orbital analysis of 
the localized bands provides an idea about the connec- 
tion between the change of ligand coordination and lig- 
and field splitting with the magnetic properties along the 
series. Fig 2(a)-(d) represent the LSDA bands of the 
relaxed structures with the Al and A2-3d localized N- 
th order muffin-tin orbital (NMTO)-downfolded bands 
[HI Qjj] (marked in blue) around Fermi-level. Localized 
3 d-band widths for these four systems are 3.71 eV, 3.76 
eV, 2.72 eV and 2.8 eV, respectively. This reduction in 
bandwidth of ~ leV from the first two to the last two 
systems can be mostly attributed to the structural dis- 
tortion. The orbital projected Co and Cu 3d-fatbands 
for systems I and IV are presented in fig. 3 and 4 respec- 
tively. The nature of orbital character distribution for 
system II and III are similar to I and IV. The labelling 
of various atoms can be seen from fig 5(a). 

Usually, for low-dimensional Cu 2+ -based dimer sys- 
tems, (2Q1I2I] eight Cu-3d a ,2_ y 2 bands localize near Er- 
in the present case, highly localized Al and A2 3d x 2 _ y 2 
bands for system I and II also implies the presence of 



dimer-interactions. For system III and IV, however, 
the dimer like localization of Al and A2- 3d x 2_ y 2 levels 
smears out due to hybridization with other orbitals indi- 
cating existence of more complex magnetic interaction. 
For Al, the bonding orbitals are x 2 — y 2 and 3z 2 — 1, 
while xy, yz and xz have both bonding and antibonding 
contributions. The situation is just reversed for A2. Al- 
3c?3 Z 2_ 1 orbitals form sigma-bond with 01 and 02 p z , 
rendering the respective bands to be the lowest in en- 
ergy for all systems. Al-3d x 2_ y 2 is pi-bonded with 01 
and 02 p x and p y and thus gets localized at a lower en- 
ergy level than xy. For Al, xz and yz orbitals are of 
highly hybridized in nature, the bonding part of which 
originates from the overlap of z-branch of xz or yz with 
01 and 02-p z . The antibonding part is from the x or 
y branch, which due to H-01 charge transfer, have little 
opportunity to interact with Ol-p x or p y . The occupied 
part of Al-3d xy originates due to the overlap with OA-p x 
and p y . For A2, 3d xy and x and y branch of 3d xz and 
2>d yz is sigma-bonded with 03-p x and p y . 3d x 2_ y 2 and 
3d 3z 2_ 1 is mostly of antibonding nature. However, for 
the third system, additional covalency with 04- p lowers 
both of the orbital energies, especially for the 3d^ z 2_ 1 . 
In the present series, competitive cationic covalency ef- 
fects are prominent from the orbital analysis. For the 
same 01 p x or p y orbitals, there is a competition be- 
tween Al- 3d x 2_ y 2, A2-3d xy and H-ls. For system I and 
II, Al-01 bond-length is less than A2-01 bond (Table 
I) indicating greater chances for Al- 3d x 2_ y 2 covalency. 
Thus for Al, 3d x 2__ y 2 bands localize at a lower energy 
compared to A2. For systems III and IV, Al-01 bond- 



FIG. 3: (colour online)The figure represents the orbitally projected NMTO fatbands for Col-d and Co2-d bands for the system 
C02PO4OH. The high-symmetry points are G, Z, T, Y, G, X, S, R, U, similar to the bands presented in fig 2. The energy-axis 
ranges from -3 to +2 eV at an interval of 0.5 eV. 




FIG. 4: (colour online)The figure represents the orbitally projected NMTO fatbands for Cul-d and Cu2-d bands for the system 
CU2PO4OH. The high-symmetry points are G, Z, T, Y, G, X, S, R, U, similar to the bands presented in fig 2. The energy-axis 
ranges from -3 to +2 eV at an interval of 0.5 eV. 



length is greater than A2-01 bond. Thus, A2- 3d x 2_ y 2 
levels stay at lower energy compared to the first two sys- 
tems. For the third system, additional lowering occurs 
due to bonding with OA-p x and p y . Also, for system III, 
strong interaction of d-bands with p x and p y orbital of 
03 and 04 and p z orbital of 01 and 02 pushed the filled 
ligand levels far below (around -4 eV) in comparison to 
the other systems(Fig 2). Investigation of hopping pa- 
rameters along relevant paths with the help of NMTO 
calculations identifies six significant couplings for such 
system, viz. the intra-tetrameric exchange J\ and J[ be- 
tween Al and A2, the dimeric exchange J2 and J' 2 along 
z-axis Ai - array and another in-plane dimeric exchange 
J3 and J a between A2. Figure 5(a) shows the exchange 
paths with arrows along with the structure and figure 
5(b) depicts the schematic diagram of predicted exchange 
couplings for this system. 



After identifying these relevant exchange paths, in the 
next step, we have estimated the exchange coupling con- 
stants and the nature of exchange along these path- 
ways by a total energy calculation for different spin- 
configurations in various supercell configurations consist- 
ing of 72 atoms. Under LSDA + U approximation, the 
total energy for each spin configuration is calculated with 
VASP-PAW potentials after ionic relaxation. The signifi- 
cant exchange coupling constants for all four systems are 
listed for the optimum U (5 eV for Cu and 4.5 eV for 
Co) and J (1 eV for Cu and 0.9 eV for Co) value in Ta- 
ble II. The insulating nature of all these four compounds 
are correctly represented with LSDA + U calculation. 
For the system I (table II), the tetrameric exchange is 
actually composed of two kinds of dimeric exchanges J\ 
and J[ along the two Col - 01 - Co2 paths (Fig 5(a)). 
Along the z-direction, the Col - 02 - Col ( J 2 ) and Col 
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FIG. 5: (colour online) (a) :The three-dimensional structure 
of one of the libethenite compound (IV) CU2PO4OH, with 
the relevant exchange paths denoted with arrows of differ- 
ent colours. The red and orange coloured arrow represents 
Ji and J((see text). The ^-direction dimeric exchanges 
are J2 (violet) and J2 (green) with in- plane dimeric exchanges 
between Cu2 J3 (magenta) and J4 (blue). Different O- 
ligands are coloured as Ol(yellow), 02(green), 03(red) and 
04(orange) (b): Schematic diagram of important exchange 
paths. 



- 01 - Col (J2) also forms a dimer. The in-plane-Co2- 
dimer comprises of Co2 - 03 - Co2 (J3) and Co2 - 01 

- 01 - Co2 (J4) exchanges. The values of Ji, J{, J2, 
J' 2 (table II) are of the same order and AFM. The Co2- 
dimeric exchange J3 is FM and J4 is AFM. Since all 
the couplings are not AFM, instead of frustration, long- 
range magnetic order prevails in this system. The order 
of magnitude of the couplings also implies that the mag- 
netic order is three dimensional which appears to be at 
par with the experimental results. 12\ With gradual sub- 
stitution of Cu in place of Col and Co2 in system II 
and III, J\ and J{ becomes more and more isotropic, 
the dimeric system will gradually turn into a tetramcric 
one. Also the ^-direction dimer turn into a monomer 
chain(J2 = J' 2 ). Although most of the exchanges for 
all the three other systems are AFM, for system III, 
exchanges along some paths are strongly ferromagnetic, 
which may lead to a long range ferromagnetic behaviour. 
For the final system, Cu 2 P0 4 OH, Cul - 01 - Cu2 ex- 
change is exactly tetrameric (Jj = J[) with 2-direction 
Cul - Cul monomer(J2 = J^)and in-plane Cu2 - Cu2 
dimers. All nearest neighbour couplings are AFM (Ta- 
ble II) as a result of geometric frustration in this partic- 



TABLE II: Relevant exchange couplings constants in meV 
for the series of mineral libethenites. 

System J 1 J[ J 2 J' 2 J3 J4 

C02PO4OH -25.84 -25.2 -22.8 -23.07 25.34 -10.56 

CuCoP0 4 OH -22.23 -24.47 -15.92 -13.6 12.43 -10.7 

Co0uPO 4 OH -18.02 -19.04 22.07 22.14 20.77 -3.62 

CU2PO4OH -20.30 -20.30 -10.12 -10.12 -9.15 -2.18 



ular system. The monomer and dimer exchange values 
are less than the tetrameric exchange implying the on- 
set of low-dimensional magnetic behaviour in this sys- 
tem. The magnetostructural details of this series are in 
resemblance with the experimental studies performed in 
reference [TT] . 

A large number of analytical investigations are avail- 
able for the spin-chain low-dimensional systems [22T 
[37] and also for more complicated systems like spin- 
trimers, tetramers.|281 12"5] octahedra [3D] etc. We may 
construct a generalized minimal Heisenberg Model spin- 
Hamiltonian describing the low-energy spin-interactions 
and excitations for the low-dimensional tetramer system 
as: 

4 n 

H = J\ Sj.Sj + J2 S2t-i-Sai+i 

n 

+ (J 3 + J4)^S 2l .S2 J+ 2 (1) 

i=l 

Here n is the number of tetramer units. For the isolated 
spin-cluster as presented in Fig 5(b), equation (1) can be 
simplified as: H — Ji(Si.S2 + Ss.Ss + S3.S4 + S4.S1) 
+ 3J 2 (S3S1) + 2J 3 (S 2 -S 4 ) + J 4 (5 2 . 5 4 ). Within first 
order approximation, the degenerate \S zto t\ — 1 triplet 
states for this system are given by 1 /\/2( \ tttl) - 1 tITT) ) , 

1/72(1 im) - 1 mi)), mt) - 1 ittt» and 

l/V2(\ im) - I mi>)- The \S zto t\ = degenerate 
singlet states are given by l/v^G llii) - I 44-tt)) and 
l/V^d flit) - I tttt))- Calculation of eigenvalues cor- 
responding to these eigenstates leads to evaluation of the 
spin-gap, i.e., the energy difference between this singlet 
ground state and triplet excited states. As a next step, 
we will investigate the effect of volume reduction on the 
magnitude of the spin-gap using this simple model Hamil- 
tonian. 



III. INVESTIGATION ON SPIN-GAP 

Among these four systems, presence of spin-gap is ex- 
perimentally confirmed only for CU2PO4OH. Hence, we 
have studied this particular system under volume reduc- 
tion (equivalent to the application of hydrostatic pres- 
sure) and size confinement. To study the volume reduc- 
tion effects, both the lattice parameters and ionic co- 
ordinates of various supercell configurations are relaxed 
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FIG. 6: (colour online) Effect of cell volume reduction is 
equivalent to application of hydrostatic pressure, upper and 
lower panel represent the comparison of lowest energy con- 
figuration for 1x1x2 and 2x1x1 supercells (details in 
text). 



with a particular magnetic arrangement. For a 1 x 1 x 2 
supercell, FM denotes the configuration where all the 
Cu-atoms are ferromagnetically alligned. In AFM1, only 
the intra-tetrameric exchange J\ is present. In AFM2, 
the monomeric exchange J 2 is also present. With grad- 
ual volume-reduction, AFM2 becomes the lowest energy 
configuration for 1 x 1 x 2 supercell, as seen at the upper 
panel of fig 6. Similarly, for a 2 x 1 x 1 supercell, AFM3 
is the magnetic configuration where Cu2 - Cu2 dimeric 
exchange J3 is present alongwith Ji . In AFM4, all of J\ , 
J3 and J4 are present. In this case, AFM4 becomes the 
lowest energy configuration(fig 6, lower panel)with vol- 
ume reduction with a very small energy difference with 
FM configuration. Thus, with gradual cell-volume re- 
duction, the monomeric and dimeric exchanges become 
more and more important in this system. In addition, we 



FIG. 7: (colour online) Upper panel depicts the variation 
of exchange couplings with cell- volume in a system where all 
four couplings are present. The lower panel represents the 
variation of spin-gap values (in K) with cell-volume as com- 
puted numerically from the simplified model Hamiltonian (2) . 
The red-line in the figure shows the equilibrium cell-volume 
and the corresponding spin-gap value ~ 150 K. 



have also observed that relaxation of the system in differ- 
ent magnetic configuration results into slightly different 
equilibrium volumes of the supercell. 

Evolution of these four effective exchange integrals 
with reduced volume, as presented at the upper panel 
of Figure 7, are computed for a 2 x 2 x 2 supercell, where 
all sorts of magnetic configurations are possible. The ex- 
change couplings are plotted with respect to an average 
cell volume of 72 atoms, to ease a comparison with Fig- 
ure 6. The equilibrium cell volume is also computed in a 
similar way. This figure indicates that with reduced vol- 
ume, all of the effective couplings become more and more 
antiferromagnetic, with J 3 and J4 orders of magnitude 
smaller than the other two. Variation of the spin-gap(in 
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TABLE III: Relevant exchange couplings constants in meV 
for the CU2PO4OH system under size-reduction. 



System 


Ji 


J'i 


■h 


J± 


J3 


Ji 


Cu 2 P0 4 OH 


-35.56 


-35.56 


20.32 


20.32 


16.407 


-10.69 



netic. So, in reduced dimension, the system may show a 
long-range ordering. 



IV. CONCLUSION 



K) with volume, as calculated with the simplified model 
spin-Hamiltonian (in text after equation 1), is also plot- 
ted at the lower panel of Fig 7. It is evident from the 
figure that with all four effective couplings present, a de- 
crease from the equilibrium cell volume (denoted by red 
dotted line in the figure) may result into an increase in 
the spin-gap magnitude. The situation is slightly differ- 
ent with an increase of cell volume, where the spin-gap 
first decreases and then increases. The marked red dot- 
ted line at the spin-gap axis indicates it's value (~ lbQK) 
corresponding to the equilibrium volume. The spin-gap 
value is quite close to the experimental observation. 1 14] 
The increase of spin-gap with volume-change is actually 
a manifestation of the variation of the couplings near 
equilibrium volume, which are shown at the inset of the 
upper panel of Figure 7. 

Effect of size-confinement is also studied after confin- 
ing a 36-atoms unit cell cluster and 72-atoms 1x1x2 
and 2x1x1 supercell clusters and calculating the effec- 
tive exchange integrals (Table III). A general observation 
not only indicates an increase of the exchange integrals 
but the geometric frustration of the bulk system seems 
to disappear as some such exchanges becomes ferromag- 



In conclusion, for the general class of Libethenite sys- 
tem, the effect of substitutional impurity induced trans- 
formation from a three dimensional antiferromagnet to 
a low-dimensional quantum spin- 1/2 tetramer system is 
investigated with the combined effects of structural and 
localized band-structure analysis along this series. The 
interplay between the structural distortion introduced by 
Cu-substitution at Co-site and its manifestation into the 
spin-magnetic structure of the series of systems are in- 
vestigated with the help of appropriate DFT techniques. 
Identification of significant exchange paths and calcula- 
tion of the corresponding exchange integrals provides the 
input to construct the model Hamiltonian for the low- 
dimensional system. All these studies have a close match 
with the earlier experimental observations. 
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